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Abstract 

In nature, microorganisms must often cope with hos- 
tile environmental conditions. To do so they have de- 
veloped sophisticated cooperative behavior and intri- 
cate communication capabilities, such as: direct cell- 
cell physical interactions via extra-membrane polymers, 
collective production of extracellular " wetting" fluid for 
movement on hard surfaces, long range chemical signal- 
ing such as quorum sensing and chemotactic (bias of 
movement according to gradient of chemical agent) sig- 
naling, collective activation and deactivation of genes 
and even exchange of genetic material. Utilizing these 



capabilities, the colonies develop complex spatio-temporal 
patterns in response to adverse growth conditions. We 
present a wealth of branching and chiral patterns formed 
during colonial development of lubricating bacteria (bac- 
teria which produce a wetting layer of fluid for their 
movement). Invoking ideas from pattern formation in 
non-living systems and using "generic" modeling we are 
able to reveal novel survival strategies which account 
for the salient features of the evolved patterns. Using 
the models, we demonstrate how communication leads 
to self-organization via cooperative behavior of the cells. 
In this regard, pattern formation in microorganisms can 
be viewed as the result of the exchange of information 
between the micro-level (the individual cells) and the 
macro- level (the colony). We mainly review known re- 
sults, but include a new model of chiral growth, which 
enables us to study the effect of chemotactic signaling on 
the chiral growth. We also introduce a measure for weak 
chirality and use this measure to compare the results of 
model simulations with experimental observations. 



1 Introduction 

Among non-equilibrium dynamical systems, living organisms 
are the most challenging ones that scientists can study. A bi- 
ological system constantly exchanges material, energy and in- 
formation with the environment as it regulates its growth and 
survival. The energy and chemical balances at the cellular level 
involve an intricate interplay between the microscopic dynamics 
and the macroscopic environment, through which life at the in- 
termediate mesoscopic scale is maintained 0. The development 
of a multicellular structure requires non-equilibrium dynamics, 
as microscopic imbalances are translated into the macroscopic 
gradients that control collective action and growth [|66| . 

Much effort is devoted to the search for basic principles of or- 
ganization (growth, communication, regulation and control) on 
the cellular and multicellular levels g g| |7|, g|, @, [73), g(| g|]. 
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Our approach is to use the successful conceptual framework for 
pattern formation in non-living systems as a tool to unravel their 
significantly more complex biological counterparts. Of critical 
importance is the choice of starting point, i.e. the choice of 
which phenomena to study: it has to be simple enough to al- 
low progress, but also well motivated by the significance of the 
results. Cooperative microbial behavior is well suited for these 
requirements, as we explain below. We focus on two examples: 
the branching growth in bacterial colonies of P. dendritiformis 
var. dendron , and the chiral growth of P. dendritiformis var. 
chiralis . Examples of these patterns are shown in Fig. |l| 

Traditionally, bacterial colonies are grown on substrates with 
a high nutrient level and intermediate agar concentration 



Such "friendly" conditions yield colonies of simple compact pat- 
terns, which fit well the contemporary view of bacterial colonies 
as a collection of independent unicellular organisms (non-interacting 
"particles"). However, bacterial colonies in nature must regu- 
larly cope with hostile environmental conditions |72, [67J. When 
hostile conditions are created in a petri dish by using a very low 
level of nutrients, a hard surface (high concentration of agar), 
or both, very complex patterns are observed. 

Drawing on the analogy with diffusive patterning in non- 
living systems |39], [47], [H| ||, U , we can state that complex pat- 



terns are expected. The cellular reproduction rate that deter- 
mines the growth rate of the colony is limited by the level of 
nutrients available for the cells. The latter is limited by the 
diffusion of nutrients towards the colony (for low nutrient sub- 
strate). Hence colony growth under certain conditions should 
be similar to diffusion limited growth in non-living systems as 
mentioned above [| [(J. The study of diffusive patterning 
in non-living systems teaches us that the diffusion field drives 
the system towards decorated (on many length scales) irregular 
fractal shapes |50|, |49|, |28j. Indeed, bacterial colonies can 

develop patterns reminiscent of those observed during growth 
in non-living systems JZg |2|, |3g, [jBJ, |T§ [j]§ |I§ g |S§ . 



But, this is certainly not the end of the story. In fact, the 
colonies exhibit far richer behavior. This is ultimately a re- 
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flection of the additional levels of complexity involved ITJ], 119 



lq, 17, 11, 12, 25, M. The building blocks of the colonies are 



themselves living systems, each having its own autonomous self- 
interest and internal degrees of freedom. Yet, efficient adapta- 
tion of the colony to adverse growth conditions requires self- 
organization on all levels which can only be achieved via co- 
operative behavior of these individual cells. Thus, pattern for- 
mation at the colony level may be viewed as the outcome of 
a dynamical interplay [14], |15|] between the micro-level (the 
individual cell) and the macro- level (the colony). For this in- 
terplay to work, the effects of changes at the micro-level must 
make themselves felt at the macro-level. This is why the notion 
of a singular perturbation, discovered to be a key for under- 
standing pattern selection in non-living systems, will be of even 
more importance here. 

This manuscript is mainly a review of our modeling studies. 
Yet we include some important new results: 1. Explanations 
and modeling of chemotaxis during chiral growth. 2. Modeling 
and characterization of weak chirality (global weak twist) during 
branching growth. 

For completeness, section includes a brief description of the 
necessary biological background needed to justify our models of 
the growth. 

How should one approach the modeling of the complex bac- 
terial patterning? With present computational power it is nat- 
ural to use computer models as a main tool in the study of com- 
plex systems. However, one must be careful not to be trapped in 
the "reminiscence syndrome", described by J. D. Cowan PS| , as 
the tendency to devise a set of rules which will mimic some as- 
pect of the observed phenomena and then, to quote J. D. Cowan 
"They say: 'Look, isn't this reminiscent of a biological or phys- 
ical phenomenon!' They jump in right away as if it's a decent 
model for the phenomenon, and usually of course it's just got 
some accidental features that make it look like something." Yet 
the reminiscence modeling approach has some indirect value. 
True, doing so does not reveal (directly) the biological func- 
tions and behavior. However, it does reflect understanding of 
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geometrical and temporal features of the patterns, which indi- 
rectly might help in revealing the underlying biological prin- 
ciples. Another extreme is the "realistic modeling" approach, 
where one constructs a model that includes in details all the 
known biological facts about the system. Such an approach sets 
a trajectory of ever including more and more details (vs. gener- 
alized features). The model keeps evolving to include so many 
details that it loses any predictive power. 

Here we try to promote another approach - the "generic 
modeling" one [||(], [18], f|, [23j. We seek to elicit, from the exper- 



imental observations and the biological knowledge, the generic 
features and basic principles needed to explain the biological 
behavior and to include these features in the model. We will 
demonstrate that such modeling, with close comparison to ex- 
perimental observations, can be used as a research tool to reveal 
new understanding of the biological systems. 

Generic modeling is not about using sophisticated, as it may, 
mathematical description to dress pre-existing understanding of 
complex biological behavior. Rather, it means a cooperative ap- 
proach, using existing biological knowledge together with math- 
ematical tools and synergetic point of view for complex systems 
to reach a new understanding (which is reflected in the con- 
structed model) of the observed complex phenomena. 

The generic models can yet be grouped into two main cate- 
gories: 1. Discrete models such as the Communicating Walkers 
models of Ben- Jacob et al. [pl|, [12], [| and the Bions model of 



Kessler and Levine |40| , |41[| . In this approach, the microorgan- 
isms (bacteria in the first model and amoebae in second) are 
represented by discrete, random walking entities (walkers and 
bions, respectively) which can consume nutrients, reproduce, 
perform random or biased movement, and produce or respond 
to chemicals. The time evolution of the chemicals is described 
by reaction-diffusion equations. 2. Continuous or reaction- 



diffusion models |)5|, |48| . In these models the microorganisms 
are represented via their 2D density, and a reaction-diffusion 
equation of this density describes their time evolution. This 
equation is coupled to the other reaction-diffusion equations for 
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the chemical fields. In the context of branching growth, this 
idea has been pursued recently by Mimura and Matsushita et 
al. |55| ], Kawasaki et al. P5| |, Kitsunezaki [H| and Ko- 
zlovsky et al. |45|. A summary and critique of this approach 



can be found in |68[ and |3T 



In section |3] we present the continuous modeling of the branch- 
ing growth. In section |] the chiral growth is modeled using 
the "atomistic" Communicating Spinors model, which enables 
us to model chemotaxis response. It is the first time that the 
chemotaxis effect on chiral growth has been studied. The ac- 
tual study of chemotaxis, both in the chiral growth and the 
branching growth is done in section |^. 

Section | is devoted to the studies of weak chirality. The phe- 
nomenon is modeled using both continuous and discrete models. 
We introduce a measure for weak chirality which enables a more 
crucial comparison between the models' results and the observed 
patterns. Good agreement was found. 

Conclusions are presented in Section [f[ We explain that 
the weak chirality phenomenon is general, and show examples 
of weak chirality during growth of the chiral morphotype and 
the vortex morphotype. In the latter it results from a different 
mechanism, and indeed the twist is not linear with the radius 
of growth. 



2 Observations and Biological back- 
ground 

Following the experimental observations which are explained in 
this manuscript, we will describe the most relevant information 
for the understanding and modeling of the observed colonial 
patterning. We base relevancy on our previous experience and 
we concentrate on bacterial movement. 



6 



2.1 Experimental observations: branching growth 
of bacterial colonies 

2.1.1 Macroscopic Observations 

Some additional examples of the patterns exhibited by colonies 
of the T morphotype are shown in figures [|, || and |j. For in- 
termediate agar concentrations (about 1.5% - 1.5g in 100m/), 
at very high peptone levels (above 10g/l) the patterns are com- 
pact (Fig. |3]a). At somewhat lower but still high peptone levels 
(about 5-10(7/7) the patterns exhibit quite pronounced radial 
symmetry and may be characterized as dense fingers (Fig. |3]b), 
each finger being much wider than the distance between fin- 
gers. For intermediate peptone levels, branching patterns with 
lower fractal dimension (reminiscent of electro-chemical deposi- 
tion) are observed (Fig. |^c). The patterns are "bushy", with 
branch width smaller than the distance between branches. As 
the peptone level is lowered, the patterns become more rami- 
fied and fractal-like. Surprisingly, at even lower peptone levels 
(below 0.25(7/7 for 2% agar concentration) the colonies revert to 
organized structures: fine branches forming a well defined global 
envelope. We characterize these patterns as fine radial branches 
(Fig. |3|d). For extremely low peptone levels (below O.lg/l), the 
colonies lose the fine radial structure and again exhibit fractal 
patterns (Fig. |J). For high agar concentration the branches are 
very thin (Fig. |]b). 

At high agar concentration and very high peptone levels the 
colonies display a structure of concentric rings (Fig. |^). At 
high agar concentrations the branches also exhibit a global twist 
with the same handedness (weak chirality), as shown in Fig. || 
Similar observations during growth of other bacterial strains 
have been reported by Matsuyama et al. |55j] . We referred 
to such growth patterns as having weak chirality, as opposed to 
the strong chirality exhibited by the C morphotype. 

A closer look at an individual branch (Fig. [F[) reveals a 
phenomenon of density variations within the branches. These 
3-dimensional structures arise from accumulation of cells in lay- 
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ers. The aggregates can form spots and ridges which are either 
scattered randomly, ordered in rows, or organized in a leaf- veins- 
like structure. The aggregates are not frozen; the cells in them 
are motile and the aggregates are dynamically maintained. 



2.1.2 Microscopic Observations 

Under the microscope, cells are seen to perform a random-walk- 
like movement in a fluid. This fluid is, we assume, excreted by 
the cells and/or drawn by the cells from the agar Jl8|, [17| . The 
cellular movement is confined to this fluid; isolated cells spotted 
on the agar surface do not move. The boundary of the fluid thus 
defines a local boundary for the branch (Fig. |8|). Whenever the 
cells are active, the boundary propagates slowly as a result of the 
cellular movement pushing the envelope forward and production 
of additional wetting fluid. Electron microscope observations 
reveal that these bacteria have flagella for swimming. 

The observations reveal also that the cells are active at the 
outer parts of the colony, while closer to the center the cells 
are stationary and some of them sporulate (form spores) (Fig. 
H). It is known that certain bacteria respond to adverse growth 
conditions by entering a spore stage until more favorable growth 
conditions return. Such spores are metabolically inert and ex- 
hibit a marked resistance to the lethal effects of heat, drying, 
freezing, deleterious chemicals, and radiation. 

At very low agar concentrations (below 0.5%) the bacteria 
swim inside the agar and not on its surface. Between 0.5% and 
1% agar concentration some of the bacteria move on the surface 
and some inside the agar. 



2.2 Chiral patterns 

Chiral asymmetry (first discovered by Louis Pasteur) exists in a 
whole range of scales, from subatomic particles through human 
beings to galaxies, and seems to have played an important role 
in the evolution of living systems [^, [| . Bacteria display var- 



ious chiral properties. Mendelson et al. |5g, |60|, Pm, [59[] showed 
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that long cells of B. subtilis can grow in helices, in which the 
cells form long strings that twist around each other. They have 
shown also that the chiral characteristics affect the structure of 
the colony. Ben- Jacob et al. |L4], [19|, [E| have found yet an- 
other chiral property - the strong chirality exhibited by the C 
morphotype. Here, the flagella handedness acts as a microscopic 
perturbation which is amplified by the diffusive instability, lead- 
ing to the observed macroscopic chirality. This appears to be 
analogous to the manner in which crystalline anisotropy leads to 
the observed symmetry of snowflakes 0] ; more about this later. 



2.2.1 A Closer Look at the Patterns 

C morphotype exhibits a wealth of different patterns accord- 



ing to the growth conditions (Fig. |l(|). As for T morphotype, 
the patterns are generally compact at high peptone levels and 
become ramified (fractal) at low peptone levels. At very high 
peptone levels and high agar concentration, C morphotype con- 
ceals its chiral nature and exhibits branching growth similar to 
that of T morphotype. 

Below 0.5% agar concentration the C morphotype exhibits 
compact growth with density variations. These patterns are 
almost indistinguishable from those developed by the T mor- 
photype. In the range of 0.4%-0.6% agar concentration the C 
morphotype exhibits its most complex patterns (Fig. [11]). Sur- 
prisingly, these patterns are composed of chiral branches of both 
left and right handedness. Microscopic observations reveal that 
part of the growth is on top of the agar surface while in other 
parts the growth is in the agar. Our model of the chiral growth 
explains that indeed growth on top of the surface and in the 
agar should lead to opposite handedness. 

Optical microscope observations indicate that during growth 
of strong chirality the cells move within a well defined envelope. 
The cells are long relative to those of T morphotype, and the 
movement appears correlated in orientation (Fig. |l^). Each 
branch tip maintains its shape, and at the same time the tips 
keep twisting with specific handedness while propagating. Elec- 
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tron microscope observations do not reveal any chiral structure 
on the cellular membrane [17|. 



2.3 Biological background 
2.3.1 Bacterial movement 

In the course of evolution, bacteria have developed ingenious 
ways of moving on surfaces. The most widely studied and per- 
haps the most sophisticated translocation mechanism used by 
bacteria is the flagellum |£7|| , but other mechanisms exist as well 



34J] . Swimming is a solitary movement done in liquid. A swim- 
ming bacterium runs nearly straight runs, interrupted by short 
periods of tumbling. Tumbling event is a random rotation in 
one location. The direction of the next run is dictated by the 
final orientation of the tumbling bacterium. 

A swimming bacterium propagates itself by rotating a bun- 
dle of flagella. Each flagellum is an helical protein filament 
which is hooked to a molecular engine transversing the bacte- 
rial membrane. The engines of all flagella rotate synchronously 
clockwise or counterclockwise. When the bacterium turn them 
counterclockwise, the flagella form an aligned bundle and push 
the bacterium forward. When they turn clockwise, the flagella 
disjoin and the bacterium tumbles. 

In our experiments, the bacteria do not swim in a pre- 
supplied liquid, but in a layer of fluid on the surface of the 
agar. The bacterial cells move individually and at random in the 
same manner as flagellated bacteria move in wet mounts (i.e., 
nearly straight runs separated by brief tumbling). Swimming 
takes place only in sufficiently thick surface fluid. Microscope 
observations reveal no organized flow-field pattern. 

Based on microscope observations of movement and electron 
microscope observations of flagella we identify the movement of 
C and T morphotype as swimming. Cells tumble about every 
tt ~ 1 — 5 sec depending on external conditions. The speed 
of the bacterium between tumbling events is very sensitive to 
conditions such as the liquid viscosity, temperature and pH level. 
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Typically, it is of the order of l-10fim/sec. 

Swimming can be approximated by a random walk with vari- 
able step size ||20|| . At low bacterial densities the random walk 
can be described by a diffusion equation with a diffusion coeffi- 
cient D b = v 2 r T = 1CT 8 — 10~ 5 cm 2 / sec. Low bacterial densities 
means that the mean free path between bacterial collisions l c is 
longer than the tumbling length It = vtt, thus collisions be- 
tween the bacteria can be neglected. The mean free path (or 
collision length) is 

in 3 dimensions ^ 
in 2 dimensions 

where p is the 3D bacterial density and o is the 2D density - 
the projection of p on the surface. 

At high densities (l c < It), the collisions cannot be neglected. 
In attempt to approximate the dynamics in those conditions, 
one may want to consider the time of straight motion to be l c j v 
instead of r T . Hence D b depends on the bacterial density to 
yield 

/ vp~* in 3D 
Db oc< _i • (2) 

[ va 2 m 2D 

This approximation is valid under the assumptions that a colli- 
sion event is identical to a tumbling event (abrupt uncorrelated 
change in direction of motion), that a tumbling event is indepen- 
dent of the collisions, and that the speed between such events 
is not affected by their frequency. 

The assumption that a collision event is like a tumbling event 
poses many problems. Even if the bacteria do not activate spe- 
cial response to collision it is unrealistic to assume that collisions 
are elastic, or that the flagella adopt immediately to the new 
orientation which changes during collisions. Thus it is reason- 
able to assume strong correlation between the cell's orientation 
before collision and the cell's orientation after collision. In ad- 
dition, the orientation after the collision should be biased ac- 
cording with the average direction of motion of the surrounding 
bacteria, as they carry the liquid with them. The important 
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parameter is not the collision length l c but re-orientation time 
T r . The re-orientation time is the time it takes a bacterium to 
loose memory of its initial orientation, i.e. the time span on 
which the final orientation has effectively no correlation with 
the initial orientation. At low densities the re-orientation time 
r r is equal to the tumbling time t t . As the density rises and the 
collisions become more frequent, r r decrease. r r defines the den- 
sities above which the constant diffusion coefficient Db = v 2 tt 
is not a good approximation. It is quite possible that these den- 
sities are high enough so as to make the velocity and even the 
type of motion dependent on bacterial density, making relation 
(0) irrelevant. In any case, high cellular densities does mean 
an effective decrease in the diffusion coefficient related to the 
bacterial movement. 

When swimming in an unstirred liquid, very low cellular 
densities also effect the movement. The bacteria secrete vari- 
ous materials into the media and some of them, e.g. enzymes 
and other polymers, significantly change the physical properties 
of the liquid making it more suitable for bacterial swimming. 
The secretion of these materials depend on cellular density, thus 
at not-too-high densities the speed of swimming rise with the 
cellular density. Hence the diffusion coefficient related to the 
bacterial movement should be a non-monotonic function of the 
bacterial density. Moreover, the specific functional form might 
depend on the specific bacterial strain. 

In other conditions there is similar but more pronounced 
effect. On semi-solid surface the bacteria cannot swim at all in- 
side the agar and they have to produce their own layer of liquid 
to swim in it. To produce such fluid the bacteria secrete lubri- 
cant (wetting agents). Other bacterial species produce known 
extracellular lubricants (such as surfactants, see p0[ [H], |26], [57] 
and references therein, or the extracellular slime produced by 
Proteus mirabilis [[f8[] ). These are various materials (various 
cyclic lipopeptides were identified) which draw water from the 
agar. The composition and properties of the lubricant of P. 
dendritiformis is not known, but we will assume that higher 
concentration of lubricant is needed to extract water from a 
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dryer agar, and that the lubricant is slowly absorbed into the 
agar (or decomposes). A single bacterium on the agar surface 
cannot produce enough fluid to swim in it, thus the bacteria 
cannot break out of the layer fluid and the branches of a T or 
C colony can be defined by this fluid. Whenever bacteria en- 
ter the shallower parts of the layer, at the edge of the branch, 
they become sluggish, indicating that the depth of the layer ef- 
fects the bacterial movement. It can be argued (see section |3~^ ) 
that in such cases the bacterial speed is related to the bacterial 
density by a power law (at least in low densities). Not only 
the diffusion coefficient related to the bacterial movement is a 
non-monotonic function of the bacterial density (as in a liquid 
agar), but it is also vanishes for extremely low densities. In this 
case it is clear that the specific functional form depend on the 
specific bacterial strain (B. subtilis, for example, cannot move 
at all under such conditions). 



2.3.2 Chemotaxis in swimming bacteria 

Chemotaxis means changes in the movement of the bacteria in 
response to a gradient of certain chemical field |], [H], H> 131 ■ 



The movement is biased along the gradient either in the gradient 
direction or in the opposite direction. Usually chemotactic re- 
sponse means a response to an externally produced field, like in 
the case of chemotaxis towards food. However, the chemotactic 
response can be also to a field produced directly or indirectly by 
the bacterial cells. We will refer to this case as chemotactic sig- 
naling. The bacteria sense the local concentration R of a chem- 
ical by membrane receptors binding the chemical's molecules 
|], |46|. It is crucial to note that when estimating gradients of 
chemicals, the bacterial cells actually measure changes in the 
receptors' occupancy and not in the concentration itself. When 
put in continuous equations [53, [JT|, this indirect measurement 



translates to measuring the gradient 

OR K OR 



dx (K + R) (K + Rf dx 



(3) 
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where K is a constant whose value depends on the receptors' 
affinity, the speed in which the bacterium processes the signal 
from the receptor, etc. This means that the chemical gradient 
times a factor Kj [K + R) 2 is measured, and it is known as the 
"receptor law" fl63| . 

In a continuous model, we incorporate the effect of chemo- 
taxis by introducing a chemotactic flux J c h em ,'- 



x(-R)V-R is the gradient sensed by the bacteria (with x(R) hav- 
ing the units of 1 over chemical's concentration). x{R) is usually 
taken to be either constant or the "receptor law", ((a) is the 
bacterial response to the sensed gradient (having the same units 
as a diffusion coefficient times the units of the bacterial density 
a). It is positive for attractive chemotaxis and negative for re- 
pulsive chemotaxis. 

Ben- Jacob et al. argued [§, [2^, [7|, |6| that for the colonial 
adaptive self-organization the bacteria employ three kinds of 
chemotactic responses, each dominant in different regime of the 
morphology diagram (the claim was made for T morphotype, 
but the same hold for their relatives C morphotype). One re- 
sponse is the food chemotaxis mentioned above. It is expected 
to be dominant for only a range of nutrient levels (see the "recep- 
tor law" below). The two other kinds of chemotactic responses 
are signaling chemotaxis. One is long-range repulsive chemo- 
taxis. The repelling chemical is secreted by starved bacteria at 
the inner parts of the colony. The second signal is a short-range 
attractant. The length scale of each signal is determined by 
the diffusion constant of the chemical agent and the rate of its 
spontaneous decomposition. 

Amplification of diffusive Instability Due to Nutrients Chemo- 
taxis: In non-living systems, more ramified patterns (lower frac- 
tal dimension) are observed for lower growth velocity. Based on 
growth velocity as function of nutrient level and based on growth 
dynamics, Ben- Jacob et al. concluded that in the case of 
bacterial colonies there is a need for mechanism that can both 




(4) 
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increase the growth velocity and maintain, or even decrease, the 
fractal dimension. They suggested food chemotaxis to be the 
required mechanism. It provides an outward drift to the cel- 
lular movements; thus, it should increase the rate of envelope 
propagation. At the same time, being a response to an exter- 
nal field it should also amplify the basic diffusion instability of 
the nutrient field. Hence, it can support faster growth velocity 
together with a ramified pattern of low fractal dimension. 

Repulsive chemotactic signaling: We focus now on the for- 
mation of the fine radial branching patterns at low nutrient 
levels. From the study of non-living systems, it is known that 
in the same manner that an external diffusion field leads to the 
diffusion instability, an internal diffusion field will stabilize the 
growth. It is natural to assume that some sort of chemotactic 
agent produces such a field. To regulate the organization of the 
branches, it must be a long-range signal. To result in radial 
branches it must be a repulsive chemical produced by bacteria 
at the inner parts of the colony. The most probable candidates 
are the bacteria entering a pre-spore stage. 

If nutrient is deficient for a long enough time, bacterial cells 
may enter a special stationary state - a state of a spore - which 
enables them to survive much longer without food. While the 
spores themselves do not emit any chemicals (as they have no 
metabolism), the pre-spores (sporulating cells) do not move and 
emit a very wide range of waste materials, some of which unique 
to the sporulating cell. These emitted chemicals might be used 
by other bacteria as a signal carrying information about the 
conditions at the location of the pre-spores. Ben- Jacob et al. 
[fL8l , |l7l , |23| suggested that such materials are repelling the bac- 
teria ('repulsive chemotactic signaling') as if they escape a dan- 
gerous location. 



2.3.3 Food Consumption, Reproduction and Starva- 
tion 

P. dendritiformis, like most bacteria, reproduce by fission of 
the cell into two daughter cells which are practically identical 
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to the mother cell. The crucial step in the cell division is the 
replication of the genetic material and its sharing between the 
daughter cells. Haste replication of DNA might lead to many 
errors - most organisms limit the rate of replication to about 
1000 bases per second. Thus the reproduction must take at 
least minimal reproduction time tr. This reproduction time tr 
is about 25mm in Bacilli. 

For reproduction, as well as for movement and other metabolic 
processes, bacteria and all other organisms need influx of energy. 
Any organism which does not get its energy directly from sun- 
light (by photo-synthesis) needs an external supply of food. In 
the patterning experiments the bacteria eat nutrient from the 
agar. As long as there is enough nutrient and no significant 
amount of toxic materials, food is consumed (for cell replication 
and internal processes) at maximal rate fl c . To estimate Q c we 
assume that a bacterium needs to consume an amount of food 
Cr of about 3 x 10~ 12 g. It is 3 times its weight - one quanta 
for doubling body mass, one quanta used for movement and all 
other metabolic processes during the reproduction time tr, and 
one quanta is for the reduced entropy of making organized cell 
out of food. Hence Q c is about 2fg/sec (lfg = 10~ 15 gram). 

If nutrient is deficient for a long enough period of time, the 
bacterial cells may enter a special stationary state - a state of 
a spore - which enables them to survive much longer without 
food. The bacterial cells employ very complex mechanisms tai- 
lored for the process of sporulation. They stop normal activity - 
like movement - and use all their internal reserves to metamor- 
phose from an active volatile cell to a sedentary durable 'seed'. 
While the spores themselves do not emit any chemicals (as they 
have no metabolism), the pre-spores (sporulating cells, see Fig. 
H) do not move and emit a very wide range of wast e materials, 
some of which unique to the sporulating cell. These emitted 
chemicals might be used by other cells as a signal carrying in- 
formation about the conditions at the location of the pre-spores. 
Ben- Jacob et al. [18 , |17| , p3[| suggested that such materials are 



repelling the bacteria ('repulsive chemotactic signaling') as if 
they escape a dangerous location. 



16 



When bacteria are grown in a petri dish, nutrients are usu- 
ally provided by adding peptone, a mixture including all the 
amino acids and sugars as source of carbon. Bacteria which are 
not defective in synthesis of any amino acid can grow also on a 
minimal agar in which a single source of carbon and no amino 
acids are provided. Such growth might seem to be easier to 
model as the growth is limited by the diffusion of a single chem- 
ical. However, during growth on minimal agar there is usually a 
higher rate of waste products accumulation, introducing other 
complications into the model. Moreover many of our strains 
are auxotrophic i.e. defective in synthesis of some amino acids 
and need an external supply of it. Providing the bacteria with 
these amino acids and only a single carbon source might pose 
us the question as to what is the limiting factor in the growth 
of the bacteria. For all those reasons we prefer to use peptone 
as nutrient source. 

We said that if there is ample supply of food, bacteria re- 
produce in a maximal rate of one division in tr. If the avail- 
able amount of food is limited, bacteria consume the maximum 
amount of food they can. In the limit of low bacterial density, 
the available amount of food over the tumbling time tt is the 
food contained in the area TTy/Df,D n , where Df, and D n are 
the diffusion coefficients of the bacteria and the food, respec- 
tively. Hence the rate of food consumption is given by ny/Di,D n 
(weather is constant or not). 

In a continuous model, reproduction of bacteria translate 
to a growth term of the bacterial density which is o times the 
eating rate per bacteria. In the limit of high nutrient it is ct/tr, 
and in the limit of low nutrient it is proportional to no. This 
brings to mind Michaelis-Menten law |k| of ^ no with K, 7 
constants. Many authors take only the low nutrient limit of this 
expression, Kno, although it is not biologically established that 
the bacteria in the experiments are limited by the availability 
of food and not by their maximal consumption rate. 
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3 Continuous models for the branch- 
ing growth of T morphotype 

3.1 The Lubricating Bacteria Model 

The Lubricating Bacteria model is a react ion- diffusion model for 
the bacterial colonies of the T morphotype |3T, This model 



includes four coupled fields. One field describes the bacterial 
density b(x,t), the second describes the height of lubrication 
layer in which the bacteria swim l(x,t), a third field describes 
the nutrients n(x, t) and the fourth field is the stationary bac- 
teria that "freeze" and begin to sporulate s(x,t) (see section 

2^3). 



We first describe the dynamics of the bacteria and of the 
nutrient. The two reaction-diffusion equations governing those 
fields are: 

— = movement + TJo, n) 
ot 

— = D n V 2 n-g(n,b) (5) 

where Tb(b,n) is the bacterial reproduction term. The nutri- 
ent diffusion is a simple diffusion process with a constant diffu- 
sion coefficient D n . The bacteria consume nutrients at the rate 
g(n, b) which is taken to be: 

g(n, b) = nb (6) 

This approximate term is correct at the limit of low nutrient 
level and low bacterial density. 

The nutrient consumed by bacteria serves as an energy source 
and as a precursor for synthesis of macromolecules. There is 
probably a minimum amount of energy necessary to maintain 
cell structure and integrity, called maintenance energy, and nu- 
trients used to supply the maintenance energy are not available 
for cell growth. We assume that those nutrients are required at 
a constant rate \x. Bacteria then cannot utilize all the nutrients 
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for reproduction: 

T b (b,n) = g(n,b) - fib (7) 

There is no explicit term for sporulation in (^). Instead the 
reproduction term T b can be used also to model sporulation. 
Sporulation is initiated by starvation so it is complementary to 
reproduction. When T b > bacteria reproduce and so sporu- 
lation is excluded. In the other case when I";, < bacteria re- 
duce in number. This should represent sporulation. Note that 
T b < when the nutrient level is low so indeed we can claim 
that bacteria are starved. This simple sporulation scheme has 
its limitations. Its rate is set by the nutrient consumption and 
effects such as density dependence are neglected. We represent 
the sporulating bacteria by the field s(x,t), whose time evolu- 
tion is: 

2l = I if r„ > o 

at \ -r 6 if T b < o l5j 

In other continuous models the bacterial reproduction and 
the sporulation process are modeled differently [§4|, [38], [53], |3lH . 
The bacterial reproduction is proportional to the nutrient con- 
sumption rate and sporulation is an independent process that 
proceeds at a rate // which could depend on other variables such 
as the bacterial density or the nutrient concentration. Note that 
if those modifications are applied to this model, the functional 
form of the bacterial reaction terms ([5]), (0) will not change. 
The equation for the field s (§) will change to: 

| = * (»> 

Since the dynamics of the other variables are separated from s, 
its different dynamics are not significant. Moreover, the modi- 
fication of the equation for s is minor. The difference is only in 
the biological interpretation of the terms. 

We now turn to the bacterial movement. In a uniform layer 
of liquid, bacterial swimming is a random walk with a variable 
step length and can be approximated by diffusion. The layer 
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of lubricant is not uniform, and its height affects the bacterial 
movement. An increase in the amount of lubricant decreases 
the friction between the bacteria and the agar surface. We sug- 
gest that the bacterial movement depends on the local lubricant 
height through a power law with the exponent 7 > 0: 

movement = V • (D b PVb) (10) 

where D b is a constant with dimensions of a diffusion coefficient. 
D b is related to the fluid's viscosity and the dryness of the agar 
might affect this viscosity. Gathering the various terms gives 
the partial model: 

V • (D h PVb) +nb-fib 
D n V 2 n - nb 

— min (nb — fib, 0) (11) 

It is possible to define dimensionless time and space variables 
t' = tfi and x' = x^J \xjD n . In those units the parameters D n 
and ji are equal to 1. 

We model the dynamics of the lubricating fluid also by a 
reaction diffusion equation. There are two reaction terms: pro- 
duction by the bacteria and absorption into the agar. The dy- 
namics of the field are: 

^ = -VJ l + f l (b,n,l)-\l (12) 

— * 

where J; is the fluid flux, fi(b, n, I) is the fluid production term 
and A is the absorption rate of the fluid into the agar. 

We assume that the fluid production depends on the bacte- 
rial density. As the production of lubricant probably demands 
substantial metabolic efforts, it should also depend on the nu- 
trient's level. We take a simple form where the production de- 
pends linearly on the concentrations of both the bacteria and 
the nutrient. The exact relation should depend on the synthetic 
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pathway of the active agents composing the lubricant. If they 
are, for example, secondary metabolites, then their production 
does not depend on the current nutrient level, but on the prior 
accumulation of primary metabolites. However, the model is not 
sensitive to the exact dependence of the lubricant production on 
the nutrient's level since the lubricant production is important 
at the front of the expanding colony, where the nutrient level is 
close to the initial level. It is reasonable that the bacteria pro- 
duce lubricant up to a height, denoted as Im, which is sufficient 
for their swimming motion. We therefore take the production 
term to be: 

f l (b,n,l) = Tbn(l M -l) (13) 

where T is the production rate. 

We turn to the flow of the lubricating fluid. The physical 
problem is very complicated. However a simplified model will 
suffice. We model the lubricant flux as a non-linear diffusion 
process: 

Ji = -DtWl (14) 

where Di is a constant with dimensions of a diffusion coefficient. 
The diffusion term of the fluid depends on the height of the fluid 
to the power v > 0. The nonlinearity causes the fluid to have 
a sharp boundary at the front of the colony, as is observed in 
bacterial colonies. The equation for the lubricant field is: 

81 

— = V • WVl) + Tbn(l M -I) -XI (15) 
at 

The functional form of the terms that we introduced are sim- 
ple and plausible, but they are not derived from basic physical 
principles. Therefore we do not have quantitative relations be- 
tween the parameters of those terms and the physical properties 
of the agar substrate. However we can suggest some relations. 
In the experiments, the agar concentration is controlled. Higher 
agar concentration gives a drier and more solid substrate. We 
shall try to find what are the effects on the lubricant layer. We 
recall that the lubricant fluid is composed from water and active 
components such as surfactants. A drier agar can increase the 
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absorption rate A. Alternatively it can diminish the amount of 
water extracted by the active components. Then either the lu- 
bricant layer will be thinner or the bacteria will have to produce 
more of the active components. The former case should decrease 
D b while the latter should decrease the production rate T. In 
both cases the composition of the lubricant fluid will change as 
the concentration of the active components will increase. The 
lubricant fluid should become more viscous, with the effect of 
Db and D\ decreasing. 

Equation (IS) together with equations ( |TTD form the Lubri- 



cating Bacteria model. For the initial conditions, we set n to 
have a uniform distribution of level no, b to be zero everywhere 
but in the center, and the other fields to be zero everywhere. 

Our results show that the model can reproduce branching 
patterns, similar to the bacterial colonies. In the experiments 
there are two control parameters: the agar concentration and 
the initial nutrient concentration. First we examine the effect 



of changing the initial nutrient concentration uq. As Fig. |13 
shows, the model produced a dense circular colony when no was 
large. The pattern became more branched and ramified as no 
decreased until no was close to 1, which is the minimal value of 
n to support growth. 

Changing the agar concentration affects the dynamics of the 
lubricant fluid. Previously we demonstrated that a higher agar 
concentration relates to a larger absorption rate A and to lower 
production rate T and lower diffusion coefficients D\ and D^. In 



Fig. [L4] we show patterns obtained with different values of the 
parameters T and A. As we expected, increasing A or decreasing 
r produced a more ramified pattern, similar to the effect of a 
higher agar concentration on the patterns of bacterial colonies. 
Similar effects are obtained by decreasing D^. 

3.2 The Non-Linear Diffusion Model 

Under certain assumptions, the Lubricating Bacteria model can 



be reduced to the non-linear diffusion model of Kitsunezaki 44 



and Cohen [22|. All the additional assumptions needed are 
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about the dynamics near zero bacterial density: 

1) The lubricant height / is much smaller than l ma x, so that the 
production of the lubricant can be assumed to be independent 
of its height 

2) The production of lubricant is proportional to the bacterial 
density to the power a > (in the simplest case taken above 
a = 1) 

3) The absorption of the lubricant is proportional to the lubri- 
cant height to the power f3 > (in the simplest case taken above 
(3=1). 

4) Over the bacterial length scale, the two above processes are 
much faster than the diffusion process, so the lubricant height 
is proportional to the bacterial density to the power of ft /a. 

5) The friction is proportional to the lubricant height to the 
power 7 < 0. 

Given these assumptions, the lubricant field can be removed 
from the dynamics and be replaced by a density dependent dif- 
fusion coefficient. This diffusion coefficient is proportional to 
the bacterial density to the power k = —2^f3/a > 0. 
The resulting model is: 

8b 

— = V(D b k Vb)+nb- fib (16) 
at 

^ = S/ 2 n-bn (17) 
dt v 1 

§ = ,* (18) 

For k > the ID model gives rise to a front "wall", with com- 
pact support (i.e. 6 = outside a finite domain). For k > 1 
this wall has an infinite slope. The model exhibits branching 
patterns for suitable parameter values and initial conditions, as 



depicted in Fig. 15. Increasing initial levels of nutrient leads to 



denser colonies, similar to the observed patterns. 
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4 The Communicating Spinors Model 
for the chiral growth of C morpho- 
type 

The Communicating Spinors Model was developed to explain 
the chirality of the C morphotype colonies. Our purpose is to 
show that the flagella handedness, while acting as a singular 
perturbation, leads to the observed chirality. It does so in the 
same manner in which crystalline anisotropy leads to the ob- 
served symmetry of snowflakes [[|. 

It is known [^7|, [79|, |74| that flagella have specific handedness. 



Ben- Jacob et al. |]T2| proposed that the latter is the origin of 
the observed chirality. In a fluid (which is the state in most 
experimental setups), as the flagella unfold, the cell tumbles 
and ends up at a new random angle relative to the original 
one. The situation changes for quasi 2D motion - motion in a 
"lubrication" layer thinner then the cellular length. We assume 
that in this case, of rotation in a plane, the tumbling has a well 
defined handedness of rotation. Such handedness requires, along 
with the chirality of the flagella, the cells' ability to distinguish 
up from down. The growth in an upside- down petri- dish shows 
the same chirality. Therefore, we think that the determination 
of up vs. down is done either via the vertical gradient of the 
nutrient concentration or via the vertical gradient of signaling 
materials inside the substrate or via the friction of the cells 
with the surface of the agar. The latter is the most probable 
alternative; soft enough agar enables the bacteria to swim below 
the surface of the agar which leads to many changes in the 
patterns, including reversing the bias of the branches. 

To cause the chirality observed on semi-solid agar, the rota- 
tion of tumbling must be, on average, less than 90° and relative 
to a specific direction. Co-alignment (orientational interaction) 
limits the rotation. We further assume that the rotation is rel- 
ative to the local mean orientation of the surrounding cells. 

To test the above idea, we included the additional assumed 
features in the Communicating Walkers model ||18||, changing 



24 



it to a 'Communicating Spinors' model (as the particles in the 
new model have an orientation and move in quasi- ID random 
walk). The Communicating Walkers model |T8| was inspired by 
the diffusion-transition scheme used to study solidification from 
supersaturated solutions |76], [77], [ftj. The former is a hybridiza- 
tion of the "continuous" and "atomistic" approaches used in the 
study of non-living systems. Ben- Jacob et al. have presented 
in the past a version of the Communicating Spinors model for 



the chiral growth The model we present here is closely 

related to a previous model of the chiral growth, but it differs 
in two features. The first is the orientation field (see below), 
which was discontinuous piecewise constant and in this model 
it is continuous piecewise linear. The second difference is the 
definition of a single run of a spinor (the stretch between two 
tumbling events), which was defined as one run per one time 
unit (i.e. each step is a run) and now is defined as variable 
number of steps in the same direction. 

The representation of bacteria as spinors allows for a close 
relation to the bacterial properties. The bacterial cells are rep- 
resented by spinors allowing a more detailed description. At the 
end of the growth of a typical experiment there are 10 8 — 10 9 
bacterial cells in the petri- dish. Thus it is impractical to incor- 
porate into the model each and every cell. Instead, each of the 
spinors represents about 10 — 1000 cells, so that we work with 
10 4 — 10 6 spinors in one numerical "experiment" . 

Each spinor has a position r*j, direction 9i (an angle) and 
a metabolic state ('internal energy') E{. The spinors perform 
an off-lattice constrained random walk on a plane within an 
envelope representing the boundary of the wetting fluid. This 
envelope is defined on the same tridiagonal lattice where the 
diffusion equations are solved. To incorporate the swimming 
of the bacteria into the model, at each time step each of the 
active spinors (motile and metabolizing, as described below) 
recalculate its direction 9[ and moves a step of size d to this 
direction. 

The direction in which each spinor moves is determined in 
two steps; first the spinor decides whether it should continue the 
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current run, that is to continue in the same direction 9[ = 9i. In 
the basic version of the model (see Sec. |5]^ for extension of the 
model) the decision is random with a specific probability p to 
continue the run. The resulting runs have geometric distribution 
of lengths, with mean run length of d/p. Once a spinor decides 
to change direction, the new direction 9[ is derived from the 
spinor's previous direction by 

9 / i = P(9 i ,®(r i ))+Ch + £ + u (19) 

Ch and £ represent the new features of rotation due to tumbling. 
Ch is a fixed part of the rotation and £ is a stochastic part, 
chosen uniformly from an interval [—i],r]} (77 constant), u is an 
orientation term that takes, with equal probabilities, one of the 
values (forward direction) or 7r (backward direction). This 
orientation term gives the spinors their name, as it make their 
re-orientation invariant to forward or backward direction. $(r*j) 
is the local mean orientation in the neighborhood of r*j. P is a 
projection function that represents the orientational interaction 
which acts on each spinor to orient 9i along the direction $(r*j). 
P is defined by 

P(a,f3) = a + (f3 - a). (20) 

Once oriented, the spinor advances a step d in the direction 
and the new location r\ is given by: 

r' { = fi + d (cos sin 9[) (21) 

The movement is confined within an envelope which is defined 
on the tridiagonal lattice. The step is not performed if r', is 
outside the envelope. Whenever this is the counter on 

the appropriate segment of the envelope is increased by one. 
When a segment counter reaches N c , the envelope advances one 
lattice step and a new lattice cell is occupied. Note that the 
spinor's direction is not reset upon hitting the envelope, thus it 
might "bang its head" against the envelope time and time again. 
The requirement of N c hits represent the colony propagation 
through wetting of unoccupied areas by the bacteria. Note that 
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N c is related to the agar dryness, as more wetting fluid must be 
produced (more "collisions" are needed) to push the envelope 
on a harder substrate. 

Next we specify the mean orientation field $. To do so, we 
assume that each lattice cell (hexagonal unit area) is assigned 
one value of $(r), representing the average orientation of the 
spinors in the local neighborhood of the center of the cell. The 
value of $ is set when a new lattice cell is first occupied by 
the advancement of the envelope, and then remains constant. 
We set the value of $(f) to be equal to the average over the 
orientations of the N c attempted steps that led to the occupation 
of the new lattice cell. The value of $ in any given point inside 
the colony is found by linear interpolation between the three 
neighboring centers of cells. Clearly, the model described above 
is a simplified picture of the bacterial movement. For example, a 
more realistic model will include an equation describing the time 
evolution of $. However, the simplified model is sufficient to 
demonstrate the formation of chiral patterns. A more elaborate 



model will be presented elsewhere [E4 



Motivated by the presence of a maximal growth rate of the 
bacteria even for optimal conditions, each spinor in the model 
consumes food at a constant rate fl c if sufficient food is avail- 
able. We represent the metabolic state of the z-th spinor by an 
'internal energy' E^. The rate of change of the internal energy 
is given by 

dEi E m 

77~ ™ J consumed j \ ) 

dt t r 
where k is a conversion factor from food to internal energy 
(K = 5 ■ 10 3 cal/g) and E m represent the total energy loss for 
all processes over the reproduction time r^, excluding energy 
loss for cell division. C consumed is C consumed = min (Q c , Q' c ) , 
where Q' c is the maximal rate of food consumption as limited 
by the locally available food (Sec. |]). When sufficient food is 
available, Ei increases until it reaches a threshold energy. Upon 
reaching this threshold, the spinor divides into two. When a 
spinor is starved for long interval of time, Ei drops to zero and 
the spinor "freezes" . This "freezing" represents entering a pre- 
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spore state (starting the process of sporulation, see section |572|). 

We represent the diffusion of nutrients by solving the diffu- 
sion equation for a single agent whose concentration is denoted 
by n(x, t): 

"777 = -^nV C bC 'consumed > (23) 
Ot 

where the last term includes the consumption of food by the 
spinors (b is their density). The equation is solved on the same 
tridiagonal lattice on which the envelope is defined. The length 
constant of the lattice must be larger than the size of the 
spinors' step d. The simulations are started with inoculum of 
spinors at the center and a uniform distribution of the nutrient. 
Both $ and the spinors at the inoculum are given uniformly 
distributed random directions. 

Results of the numerical simulations of the model are shown 
in Fig. pl| These results do capture some important features of 
the observed patterns: the microscopic twist Ch leads to a chiral 
morphology on the macroscopic level. The growth is via stable 
tips, all of which twist with the same handedness and emit side- 
branches. The dynamics of the side-branches emission in the 
time evolution of the model is similar to the observed dynamics. 

For large noise strength 77 the chiral nature of the pattern 
gives way to a branching pattern (Fig. [I7|). This provides a 
plausible explanation for the branching patterns produced by 
C morphotype grown on high peptone levels, as the cells are 
shorter when grown on a rich substrate. The orientation inter- 
action is weaker for shorter cells, hence the noise is stronger. 



5 The effect of chemotaxis 

So far, we saw that the models can reproduce many aspects 
of the microscopic dynamics and the patterns in some range of 
nutrient level and agar concentration, but at least for the T 



-like growth, other models can do the same |31|, and reference 
there in]. We will now extend the Non-Linear Diffusion model 
and the Communicating Spinors model to test for their success 
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in describing other aspects of the bacterial colonies involving 
chemotaxis and chemotactic signaling (which are believed to by 
used by the bacteria [§, 0, |§])- 



5.1 Chemotaxis in the Non-Linear Diffusion 
Model 



As we mentioned in section p.3.2j , in a continuous model we in- 



corporate the effect of chemotaxis by introducing a chemotactic 

flux J c hem- 

Jche m = ((b)x(R)VR (24) 

x(i?)V-R is the gradient sensed by the bacteria (with x(R) hav- 
ing the units of 1 over chemical's concentration). x{R) is usu_ 
ally taken to be either constant or the "receptor law", ((b) is 
the bacterial response to the sensed gradient (having the same 
units as a diffusion coefficient times the units of the bacterial 
density b). In the Non-Linear Diffusion model the bacterial dif- 
fusion is Db = D b k , and the bacterial response to chemotaxis 
is ((b) = ( b (^D b k ^ = ( D b k+1 . ( is a constant, positive for 
attractive chemotaxis and negative for repulsive chemotaxis. 

We claimed that the fine radial branching patterns at low 
nutrient levels result from repulsive chemotactic signaling. The 
equation describing the dynamics of the chemorepellent contains 
terms for diffusion, production by pre-spores, decomposition by 
active bacteria and spontaneous decomposition: 

— = D R V 2 R + sT R - n R bR - X R R (25) 

where D R is a diffusion coefficient of the chemorepellent, T R is 
an emission rate of repellent by pre-spores, fl R is a decomposi- 
tion rate of the repellent by active bacteria, and X R is a rate of 
self decomposition of the repellent. 

Fig. [0| demonstrates the effect of repulsive chemotactic sig- 
naling. In the presence of repulsive chemotaxis the pattern be- 
comes much denser with a smooth circular envelope, while the 
branches are thinner and radially oriented. 
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5.2 Chemotaxis in the Communicating Spinors 
Model 



The colonial patterns of C morphotype (e.g. Fig. [16]) are 
rarely as ordered as the simulated patterns of the Communicat- 
ing Spinors model. For example, the branches of the observed 
colonies usually have varying curvature. In the simulations of 
C morphotype shown in Fig. |TB] all the branches have a uni- 
form curvature. One of the reasons for this difference is the 
simplifications taken during the model's development. A more 



elaborate model that we will present p4j will be a better de- 



scription of the colony. However, some of the observed features 
can be explained in the context of the Communicating Spinors 
model. In some of the observed patterns (Fig. [TTJb) , the curva- 
ture of the branches has a distinct relation to the branch's radial 
orientation (the orientation relative to the radial direction): the 
curvature is smaller when the branch is in the radial orientation 
and larger when the branch is orthogonal to that orientation. 
This brings to mind the radial organization of branches in the 
T morphotype, and indeed we where able to explain the chiral 
pattern with the aid of the same concept - repulsive chemotaxis. 

Chemotaxis was introduced in previous versions of the Com- 
municating Walkers model by varying, according to the chemi- 
cal's gradient, the probability of moving in different directions 
[ |18| , |J. Modulating the directional probability is not the way 
bacteria implement chemotaxis - they modulate the length of 
runs. However, the growth of T morphotype is insensitive to 
the details of the movement. Modulating the directional proba- 
bility is as good an implementation of chemotaxis as many other 
implementations (it was chosen for computational convenience). 
The pattern of the C morphotype is based on amplification of 
microscopic effects (singular perturbation) such as the left bias 
in the bacterial tumbling. Small differences in the microscopic 
dynamics of chemotaxis might affect the global pattern. Indeed 
we found that modulating the directional probability yield un- 
realistic results in the simulations of C morphotype. We had to 
resort to the bacterial implementation of chemotaxis - modu- 
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lating the length of runs according to the chemical's gradients. 

When modulating the length of runs of walkers or spinors 
one must be careful not to change the particles' speed. Such 
change is not observed in experiments [[H], [H] and it has far 
reaching effects on the dynamics. Changing the particles' speed 
is like changing the diffusion coefficient of the bacterial density 
field, a change that can have undesirable effects on the pattern. 

Modulating the length of spinors' runs without changing 
their speed can be done by modulating the number of steps 
that compose a single run (that was our motivation for divid- 
ing the runs into steps). Since the mean number of steps in a 
run is determined by the reorientation probability p, chemotaxis 
should modulate this probability. For chemotaxis, the probabil- 
ity of changing direction by the i-spinor in one time step is (for 
a repellent R): 

p*=p + X (R)d ei R (26) 

where R is measured at the spinor position r*j, x{R) is ^ ne same 
as in the continuous model (either constant or the "receptor 
law") and dg i is directional derivative in the spinor's direction 
0j. p* is truncated to within the range [0, 1] as it is a probabil- 
ity. The length of the resulting runs will depend on the runs' 
direction, where a spinor moving up the gradient of the repellent 
will have shorter mean run length than the same spinor moving 
down the gradient. 

The production and dynamics of the repulsive chemotactic 
signaling in the Communicating Spinors model is the same as 
in the Non-Linear Diffusion model, see Eq. fl25|) (with s rep- 
resenting the density of spinors that "freezed"). The patterns 
resulting from including repulsive chemotaxis in the model have 
indeed branches with variable curvature, as can be seen in Fig. 
[Tj|. The curvature is smaller for branches in the radial direc- 
tion. Food chemotaxis also varies the branches' curvature, but 
in a less ordered manner, not similar to the observed bacterial 
patterns. 

Under different growth conditions the C morphotype can 
produce very different patterns. As mention above, if the agar 
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is soft enough the bacteria can move inside it . In such case, the 
bias in the bacterial movement might change or even reverse, 
and it is manifested in the curvature of the branches. Widely 
changing curvature of the branches can be seen in Fig. |TT|a. 
The agar hardness was tuned such that in the beginning of the 
growth the bacteria could swim inside the agar, but they are 
forced to swim on the agar by the end of the growth due to 
the marginal water evaporation during the growth. In Fig. |20| 
we demonstrate the models' ability to explain such patterns by 
changing the spinors' bias Ch during the simulation. Ch is set 
to be a continuous random function of the colonial size, which 
is constrained only at the beginning and end of growth to have 
certain values. The function for Ch is the same in all the im- 
ages of Figs. [Zl], only that in various types of chemotaxis are 
used. As can be seen, repulsive chemotactic signaling is needed 
to explain the observed bacterial patterns. 



6 Weak chirality in T morphotype 

Colonies of T morphotype grown on hard substrate (above 2.0% 
agar concentration) exhibit branching patterns with a global 
twist with the same handedness, as shown in Figs. ^ and ^1 . 
Similar observations during growth of other bacterial strains 
have been reported by Matsuyama et al. |56|, ^5[. We refer to 
such growth patterns as having weak chirality, as opposed to 
the strong chirality exhibited by the C morphotype. 

In [JT3], Ben- Jacob et al. proposed that, in the case of T mor- 
photype, it is the high viscosity of the "lubrication" fluid during 
growth on a hard surface that replaces the cell-cell co-alignment 
of the C morphotype that limit the rotation of tumbling. They 
further assumed that the rotation should be relative to a speci- 
fied direction. They used gradient of a chemotaxis signaling field 
(specifically, the long-range repellent chemotaxis) as a specific 
direction, rather than the local mean orientation field which is 



used in the case of C morphotype. It was shown in [0 that 



inclusion of the above features in the Communicating Walkers 
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model indeed leads to a weak chirality which is highly remi- 
niscent of the observed one. The idea above also provides a 
plausible explanation to the observations of weak chirality by 
Matsuyama et al. |56| in strains defective in production of "lu- 
brication" fluid. 

6.1 Weak chirality — the Non-Linear Diffu- 
sion Model 

In the reaction-diffusion model, weak chirality can obtained by 
modifying the chemotactic mechanism and causing it to twist: 
We alter the expression for the chemotactic flux J c hem (Eq. £|) so 
that it is not oriented with the chemical gradient (VR) anymore. 
Instead it is oriented with a rotated vector R(#) Vi2, where R(#) 
is the two-dimensional rotation operator and 9 is the rotation 
angle. The chemotactic flux is thus written: 

J ch e m = ((b)x(R)(R(0)VR) (27) 

The effect of rotating the repulsive chemotaxis, as depicted in 



Fig. |22], is to make the pattern chiral, with the degree of chirality 
determined by the rotation angle 9. 

One must note that adding a similar rotation to the food 
chemotaxis does not have the same effect, because the nutri- 
ent gradients do not posses the long-range, radial nature of the 
chemorepellent gradients. 

6.2 Weak chirality — the Communicating Spinors 
Model 

As was demonstrated in section |], the Communicating Spinors 
model is robust enough to reproduce patterns of T morphotype, 
as well as patterns of C morphotype. Here we use it to model 
patterns of weak chirality. Two simulated T -like colonies are 
shown in Fig. |23l. Fig. E3k shows colony with radial branches 



while Fig. [23| b shows colony with weak chirality and thinner 



branches. In the two simulations the spinors have exactly the 
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same response to chemotaxis and the same bias Ch = 9°. The 
two runs differ in the freedom of rotation 77; in 23a the spinors 
can rotate freely (rj = 180°) while in |2"B| b the spinors rotation 
is somewhat limited (rj = 35°, while for the colony of strong 
chirality rj = 5°). 

It seems that both models - the Non-Linear Diffusion model 
and the Communicating Spinors model - can capture the es- 
sential features of the observed weak chirality Yet a closer 
examination reveals that the description of the two models is 
incompatible. In the Non-Linear Diffusion model the bias from 
the direction of the gradient is through the chemotaxis process. 
The spinors, like the bacteria, cannot modulate their runs as 
a function of the difference between their direction and that of 
the gradient; they do not know what is the direction of the gra- 
dient, only the directional derivative along their path. As was 
demonstrated in figure |23|, one of the key features for the weak 
chirality in the spinors model is the correlation in orientation of 
the spinors (through $). In fact in this model the twist of the 
branches stems from the deviation of the runs' direction from 
the orientation of neighboring spinors. The twist of the branches 
is related only indirectly, through the neighbors' orientation, to 
the chemorepellent's gradient. A continuous model of such pro- 
cess should include information about the mean orientation of 
the bacterial cells. It should include chemotaxis without rota- 
tion, anisotropic diffusion (smaller diffusion coefficient in orien- 
tations orthogonal to the mean orientation of neighbors) and a 
rotation on the diffusion operator. Such a model will be pre- 
sented elsewhere p4| . 

The discrete spinors model allows for a detailed represen- 
tation of the bacterial properties. The macroscopic dynamics 
and resulting patterns, however, are similar in both models - 
apparently the growth does not amplify the difference in the 
microscopic dynamics. Thus the 'unrealistic' microscopic de- 
scription of the Non-Linear Diffusion model does not rule it 
out as an approximation to the growth dynamics of bacterial 
colonies with weak chirality. 
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6.3 Chirality measure 



All through this manuscript we referred to chirality as a two- 
valued property - either the pattern is chiral or it is not. There 
are various attempts in the literature to quantify chirality with 
a continuous measure. See for example the method of Avnir et 
al. P7[ , who applied their method also to large disordered 
objects. While this method is general and can quantify with a 
single number the measure of asymmetry of any given object, 
it is not quite satisfactory for our purpose. We would like to 
know the time evolution of the chirality of a colony, and not 
just "mean" chirality given by a single number. We sacrifice 
the generality of the measure to that end. 

Since the growth velocity of the colonies (both experimental 
and simulated) is constant, we measure the chirality as a func- 
tion of radius instead of function of time. Thus we can work 
on chirality of an image, not of a process. The image can be 
a scanned picture of the real colony or a result of a computer 
simulation. We look for a mapping of the image to a new one, 
which in some sense does not distinguish left from right (the 
ambiguity stems from the fact that a large random object will 
not have, in general, reflection symmetry, thus there is no trivial 
definition for chirality of such objects). The mapping is defined 



where each point in the image is described by the polar coor- 
dinates (r, 8), measured from the center of the colony. Thus, 
each point is rotated by an increment A8 which depends on the 
radius r (i.e. the distance from the center). 

Working on many experimental patterns, as well as sim- 
ulated patterns, we have learned that in most cases a linear 
dependence of A9 on r is sufficient to give quite satisfactory 
results, that is, to transform a chiral pattern into a "normal" 
branching pattern. The rotating angle is thus written: 



by: 






max 



(29) 
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where r max is the radius of the colony, and 6 max is the rotation 
angle at that radius. 

The fact that this linear angular mapping suffices to "de- 
chiral" the simulated patterns may not be of much importance 
(in the case of the continuous model, at least, this is almost a 
direct result of the way in which we introduce the weak chiral- 
ity). The same transformation works for images of real colonies 
of T morphotype, but does not work for chiral colonies of other 
bacteria (see Sec. 0). This strengthens our belief in the models. 

7 Conclusions 

We first briefly reviewed experimental observations of colonial 
patterns formed by bacteria of the species Paenibacillus dendri- 
tiformis. We described both the tip-splitting growth of the T 
morphotype and the chiral, twisted-branches growth of C mor- 
photype . Both colonial patterns and optical microscope obser- 
vations of the bacteria dynamics were presented. 

In this manuscript we presented observations of various forms 
of chiral patterns in bacterial colonies. Our goal was to explain 
the various aspects of chirality. We used two types of models: 
continuous reaction-diffusion models which deal with bacterial 
density, and a hybrid semi-discrete model which deals with prop- 
erties of the individual bacterium. From a comparison of the 
models' simulation and experimental observations we conclude 
that chemotactic signaling plays an important part in the de- 
velopment of colonies of the two types. We also estimate how 
sensitive the growth is to the details of the microscopic dynam- 
ics, demonstrating that the more 'complex' the pattern is, the 
more sensitive the growth is to the small details. 

We would like to note that the P. dendritiformis is not the 
only bacteria whose colonies exhibit chirality. Ben- Jacob et al. 
discussed in |J the formation of colonies of Paenibacillus vortex, 
where each branch is produced by a leading droplet and emits 
side branches, each with its own leading droplet. Each leading 
droplet consists of hundreds to millions of bacterial cells that 
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circle a common center (a vortex) at a cellular speed of about 
10/im/s. In Fig. 24 we show a colonial pattern of these bacteria. 
The chirality we termed 'weak chirality' is evident in this figure. 
In this case the chirality is not related to the handedness of the 
flagella, but to the rotation of the vortices. When "pushed" 
by repulsive chemotaxis, Magnus force acts on the vortices and 
drive them side-ways from the radial direction of the chemore- 
pellent's gradient. This difference in mechanisms is expressed in 
the global pattern: the colonial patterns of P. vortex cannot be 
"de-chiraled" by the transformation (linear angular mapping) 
that "de-chiral" the T morphotype. The fact that the models 
for weak chirality match in this respect the the weak chiral- 
ity of T morphotype and not the 'weak chirality' of P. vortex 
is another support for their success in describing the bacterial 
colonies. 

We hope we have convinced the reader that chirality in pat- 
terns of bacterial colonies gives important clues about the un- 
derlying dynamics. The processes leading to such patterns are 
more complex than those leading to non-twisted branching pat- 
terns. The chiral patterns are more sensitive to the underlying 
dynamics and as such they require more accurate models. This 
reflects on the success of the models we presented as being a 
good description of the colonies. 
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Figure Captions 

1. (a) Typical example of branching growth of the Paeni- 
bacillus dendritiformis var. dendron (referred to as T mor- 
photype |IJ)for 0.5 g/1 peptone and 1.75 % agar concen- 



tration, (b) Chiral growth of the Paenibacillus dendriti- 
formis var. chiralis (referred to as C morphotype |IIJ) for 
2.5 g/1 peptone and 1.25 % agar concentration. 

Patterns exhibited by the T morphotype as function of 
peptone level (increasing from left to right) and agar con- 
centration (1.5% bottom row, 2% middle row, 2.5% top 
row). 

Examples of typical patterns of T morphotype for inter- 
mediate agar concentration, (a) Compact growth for 12g/l 
peptone level and 1.75% agar concentration, (b) Dense 
fingers for 3g/l peptone and 2% agar, (c) Branching frac- 
tal pattern for lg/l peptone and 1.75% agar, (d) A pattern 
of fine radial branches for O.lg/l peptone and 1.75% agar. 

Colonial patterns of T morphotype . (a) Fractal pattern 
for O.Olg/l peptone level and 1.75% agar concentration, 
(b) Dense branching pattern for 4g/l peptone and 2.5% 
agar. Note that the branches are much thinner than those 
in Fig. ^b, i.e. the branches are thinner for higher higher 
agar concentrations. 

(a) Pattern of concentric rings superimposed on a branched 
colony for 2.5g/l peptone level and 2.5% agar concentra- 
tion, (b) Concentric rings in a compact growth for 15g/l 
peptone level and 2.25% agar concentration. 
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6. Weak chirality exhibited by the T morphotype during 
growth on 4g/l peptone and 2.5% agar concentration. 

7. Density variations within branches of a colony of T mor- 
photype . Optical microscope observations x50. 

8. Closer look on branches of a colony, a) x20 magnification 
shows the sharp boundaries of the branches. The width of 
the boundary is in the order of micron, b) Numarsky (po- 
larized light) microscopy shows the hight of the branches 
and their envelope. What is actually seen is the layer of 
lubrication fluid, not the bacteria, c) x50 magnification 
shows the bacteria inside a branch. Each bar is a single 
bacterium. There are no bacteria outside the branch. 

9. Electron microscope observation of T bacteria. Round 
or oval shapes with bright center are spores. Elongated 
shapes are living cells. The cells engulfing oval shapes are 
pre-spores. 

10. Patterns exhibited by the C morphotype for different growth 
conditions, a) Thin disordered twisted branches at 0.5g/l 
peptone level and 1.5% agar concentration, b) Thin branches, 
all twisted with the same handedness, at 2g/l peptone 
level and 1.25% agar concentration, c) Pattern similar to 
(b) but on softer agar: 1.4 g/l peptone level and 0.75% 
agar concentration, d) Four inocula on the same plate, 
conditions of lg/l peptone level and 1.25% agar concen- 
tration. 

11. In agar soft enough for the bacteria to swim in it, the 
branches lose the one-side handedness they have on harder 
agar. The two colonies of (a) and (b) are of 5g/l peptone 
level and 0.6% agar concentration. The two patterns are 
of two stains of the C morphotype, strains whose patterns 
are indistinguishable on harder agar, c) Closer look (xlO 
magnification) on a colony grown at 8g/l peptone level 
and 0.6% agar concentration. 
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12. Optical microscope observations of branches of C morpho- 
type colony, a) x20 magnification of a colony at l.Qg/l 
peptone level and 0.75% agar concentration, the ant i- clockwise 
twist of the thin branches is apparent. The curvature of 
the branches is almost constant throughout the growth. 

b) xlO magnification of a colony at Ag/l peptone level 
and 0.6% agar concentration at branches are not thin, but 
have a feathery structure. The curvature of the branches 
varies, but it seems that at any given stage of growth the 
curvature is similar in all branches. That is, the curvature 
is a function of colonial growth, c) x500 magnification of 
a colony at 1.6g/l peptone level and 0.75% agar concen- 
tration. Each line is a bacterium, the bacteria are long 
(5-50/im) and mostly ordered. 

13. Effect of varying the initial nutrient concentration no on 
colony pattern. n increases from left to right: 1.2 (a), 1.4 
(b), 1.7 (c), 2 (d), 3 (e), 6 (f). The minimal value of n to 
support growth is 1. 

14. Effect of varying A, the fluid absorption rate, on colony 
pattern. The fluid production rate T is 1 in the upper 
row and 0.3 in the lower row. In both rows A increases 
from left to right: A = 0.03 (left), A = 0.1 (center), A = 1 
(right) The patterns become more ramified as A increases. 
Decreasing T also produces a more ramified pattern. The 
other parameters are: D b = D\ = l,n = 1.5 

15. Growth patterns of the Non-Linear Diffusion model, for 
different values of initial nutrient level hq. Parameters are: 
-Do — 0.1, k — l,u = 0.15. The apparent 6-fold symmetry 
is due to the underlying tridiagonal lattice. 

16. A morphology diagram of the Communicating Spinors model 
for various values of iV c and initial n concentration n . 
Ch = Q°, r? = 3°, d = 0.2, p = 0.5. 

17. When the noise rj is increases to rj = 180° the tumbling 
of the spinors becomes unrestricted. Their movement be- 
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comes like that of the T bacteria and accordingly the sim- 
ulated colonial pattern is like that of T morphotype. On 
the left 77 = 3°, on the right rj = 180°. 

18. Growth patterns of the Non-Linear Diffusion model with 
food chemotaxis (left, see section 0) and repulsive chemo- 
tactic signaling (right) included. Xof — 3,Xoi? — LAr — 
1,Tr = 0.25, Qr = 0,Ar = 0.001. Other parameters are 
the same as in figure [0^. The apparent 6-fold symmetry 
is due to the underlying tridiagonal lattice. 

19. The effect of repulsive chemotactic signaling on the Com- 
municating Spinors model, a) Without chemotaxis. b) 
With repulsive chemotaxis. The Spinors are repelled from 
the inner parts of the colony. The resulting curvature 
of the branches is reduces when they are in the radial 
direction. In spite the reversed handedness, the pattern 
resemble Fig. |TT|b. 



20. The snake-like branches observed in Fig. 11a can be re- 
produced by the Communicating Spinors model. Ch is a 
continuous function of the colony's radius (the same func- 
tion in a, b, and c). Maximal value of Ch is 8°, minimal 
value is —2°. (a) With repulsive chemotactic signaling, 
(b) Without chemotaxis. (c) With food chemotaxis. The 
best resemblance to the observed colony is obtained with 
repulsive chemotactic signaling. 

21. Weak chirality (global twist of the branches) exhibited by 
the T morphotypefor a peptone level of 0.25(7 / Z peptone 
level and agar concentration of 1.75 %. 

22. Growth patterns of the Non-Linear Diffusion model with 
a "squinting" repulsive chemotactic signaling, leading to 
weak chirality. Parameters are as in the previous picture, 

e = 43°. 

23. Weak chirality of the T morphotype is modeled by the 
Communicating Spinors model. Both simulations are with 



49 



repulsive chemotactic signaling and with bias in the walk- 
ers rotation, a) With free rotation (rj = 180°) The pat- 
tern is branched, without apparent chirality. expressed, 
b) With constrained rotation (rj = 35°) weak chirality is 
expressed. 

24. A colony of Paenibacillus vortex on 10 g/l peptone level 
and 2% agar concentration. The dots at the tips of the 
branches are bacterial vortices - each is composed of up 
to millions of bacterial cells rotating around a common 
center. The twist of the branches results from a Magnus 
force induced by repulsive chemotactic signaling. 
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